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\^ ' Using a path integral approach we rederive a recently found representation 

■ of the Casimir energy for a sphere and a cylinder in front of a plane and derive 

the first correction to the proximity force theorem. 

Oh! 

^ : 1 Introduction 



In view of the rapid experimental progress in the Casimir force measurements and in- 
creasing interest in the field of nanotechnology there is a growing demand for reliable 
I theoretical predictions. Beyond geometries allowing for a separation of variables one 

was left so far basically with the proximity force approximation for ideal reflecting 
surfaces or with the pairwise summation for dilute dielectric bodies. Also the basic 
formula for the interaction of dielectric bodies, the Lifshitz formula, is known for sim- 
ple geometries only (for a review see [1]). There is a number of approximate methods 
trying to overcome these limitations. The oldest one is the semiclassical approach 
(see, for example the recent papers [2], [3] and references therein) which exists in 
several modifications including the optical path approximation |4|. Another method 
is the multiple reflection approximation which recently has been reconsidered in the 
context of the Casimir effect [3]. New developments are the multi-scattering ap- 
proach of [6] which was applied to the geometry of two spheres and of a sphere in 
front of a plane and the calculation in [7] for a cylinder and a plane using a path 
integral formulation together with a trace formula. Significant progress had been 
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reached in the world hne method [8] which is now able to provide corrections to the 
small separation behavior. 

In the present paper we use the path integral approach developed in [9] and 
[To] and apply it to the geometry of a sphere and a cylinder in front of a plane for 
Dirichlet boundary conditions. We rederive the formulas obtained in [6] and [7j in 
a simpler way and calculate the first correction to the proximity force theorem. 

The paper is organized as follows. In the next section we describe the method of 
[9] and [10] in detail. In the third section we apply the method to the geometry of 
a sphere in front of a plane. In the fourth section we discuss some other geometries. 
Conclusions are given in the final section. Throughout the paper we put h = c = 1 



2 Functional integration and boundary conditions 

A simple, but powerful method to introduce boundary conditions into a quantum 
field theory is to use functional delta functions. This method was introduced in 
[9| for the calculation of the radiative corrections to the Casimir effect for parallel 
plates in a covariant gauge. We present it here in the simpler case for a scalar 
field theory with Dirichlet boundary conditions. The basic idea of the method is to 
restrict the space of fields over which the functional integration runs to such that 
fulfil the boundary conditions. Technically this can be achieved by representing the 
generating functional of the Green's functions as 



Z{J) = j Dcj) W5{(l){x)) exp{-S} 



with an action 



S= I dx ( ^(p{x) (-(9^ + m^) (j){x) + J{x)(j){x) ) , (2) 



where we included the term with the source J{x). This is the usual action for a 
scalar field with mass m. For simplicity we write all formulas in the Euclidean 
formulation. 

The functional Z{J), ([T]), taken without the delta function is the generating func- 
tional of the corresponding field theory without boundary conditions. Obviously, in 
([1]) more complicated actions can be inserted, especially such which include inter- 
actions. However, in connection with the Casimir effect we restrict ourselves here 
to a free field theory, i.e., to a quadratic action. Also we mention that the initial 
space where the field is defined does not need to be fiat or free of boundaries. 
The method to implement boundary conditions through functional delta functions 
is completely general, at last on the formal level taken in this section. 

The next step in the method is the use of a Fourier representation for the delta 
function, 

llS{(P{x)) = I Dh exp I dz h{z)<PU{z))\ , (3) 
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where we described the surface by a function f{z) such that x = f{z) E S and 
the coordinates z are a suitable parameterization of S. In fact, ([3]) is a functional 
integral over a scalar field h{z) defined on the surface 5*. Using ([3]) in the generating 
functional ([I]) we represent Z{J) in the form 



Z{J) = J D(j) Db exp(-5) (4) 

with 

S = ^jdx (-9^ + m^) 

—i dx / dz (f){x)H{x^ z)b{z) + / dxJ{x)(j){x) . (5) 



J Js J 
Here we introduced the new notation, 

H{x,z) = 6{x-f{z)) (6) 

projecting the argument x on the surface 5*. Obviously, S is quadratic in the fields 
(f){x) and b{z), hence S can be diagonalized, 

S = ^jdx - M^)) i-d' + m') - M^)) 

+ \l dz I dz' {b{z)-bo{z))K{z,z'){b{z)-bo{z)) 
1 /■ , /• . . s 



with 



- I dx I dyJ{x) 'D{x,y)J{y) (7) 



Mx) = J dx' D{x,x') l^t J dz H{x',z) b{z) - J{x') j , 
bo{z) = -i J dz' j dxdx' K-\z,z')H{x,z')D{x,x')J{x'), 



(8) 



and 



K{z,z') = J dx J dy H{x, z)D(x,y)H{y, z') 

= D{f{x),f{y)). (9) 

The function D{x,y) is the inverse of the kernel of the free action, 

{-d^ + m^) D{x,y) = S{x~y). (10) 
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It is just the usual propagator of the free field theory defined by the action S, ([2]). 
Instead, the new quantity 

'D{x,y) = D{x,y)~ [ dz [ dz' D{x, f{z))K~\z, z')D{f{z'),y), (11) 

Js Js 

which we introduced in ([7]), is the propagator of the given field theory with boundary 
conditions. Finally, K~^{z, z') is the inverse of K{z, z') on the surface S, 



I dz' K{z, z')K~\z\ z") = 6{z - z") 

J s 



12) 



This inversion and further operations like the trace or the determinant with K{z, z') 
are to be understood as done with an operator K whose integral kernel in a space 
of suitable functions defined on the surface S is just K{z, z'). 

With Eq.( fT2|) it is easy to check that ^D{x,y), f fTTl) . is indeed the propagator 
which fulfills the boundary conditions. First we remark that it is a propagator. 
Obviously equation (fTO!) is fulfilled for x ^ S. Inserted into Eq. lfTOj) . the first term 
delivers the necessary delta function. The second term delivers a delta function too 
which is, however, non vanishing for x & S only. This is not a problem because 
in the presence of boundary conditions the equation has to be fulfilled outside the 
boundary only. So it remains to check that the boundary conditions are fulfilled. 
Put X on the surface S, i.e., consider x E S. Then in the second term we can apply 
f|T2|) and we get zero as required by the Dirichlet condition. 

The calculation of the inverse on 5 of K is the main step in the method. All 
steps before are merely a formal rewriting. The inversion of K on S, however, is 
equivalent to solve the wave equation with boundary conditions and, in this sense, 
not much progress is archived. 

After rewriting the exponential in the form given by Eq. (j7]), it is quadratic in 
the fields. The integrals are Gaussian and after integration we represent Z{J), ([Tj), 
in the form 

Z{J) = C (det (9^ +m2))~^ (detir(2,/))"^ 

xexpj-^ / / ^^^^^) '^/^(x,y)J(y)|, (13) 

where C is a constant. In this formula the first determinant comes from the in- 
tegration over 0. It is the usual one appearing also without boundary condition. 
Of course, it is independent on the boundary conditions. The second determinant 
comes from the integration over h. It does depend on the boundary conditions. The 
source term has the standard form and it generates the usual Feynman rules with 
the boundary dependent propagator ^D{x,y) instead of the usual one. In this way 
the field theory with boundary conditions was constructed in [9]. 

The last step in the method, which was done in |^10j, is to consider the free energy 

F = -llnZ(0), (14) 
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where T is the total time which is the logarithm of the generating functional at 
zero source term. In [lOj this was done for a theory with finite temperature, but of 
course we can use it at zero temperature too. In this way we get for the boundary 
dependent part of the Casimir energy 

^casimir = - ^ In {det K {z , z'))' = ^ trlnir. (15) 

The trace in this formula is to be taken in the space of functions on the boundary 
surface S and K is the operator in this space, whose integral kernel is the function 
K{z,z'), defined in Eq.Q. In this way the Casimir energy is represented by an 
expression on the surface S. The contribution of the free space, which contains the 
highest ultraviolet divergences, had been dropped in ( ITSi) . 

Let us mention that the evaluation of the trace in ( fT5|) is a problem of the same 
difficulty as the inversion of K or the solution of the wave equation with boundary 
condition, and again, so far we have only a formal rewriting of the original problem. 



3 A sphere and a cylinder in front of a plane 

In this section we apply the general formula (fT5|) for the Casimir energy to the 
geometry of a sphere and of a cylinder in front of a plane. We start with the remark 
that the derivation in the preceding section is valid also if the initial space is not 
the full but only a half of it with z < a and if the field 0(x) satisfies Dirichlet 
(or other) boundary conditions on the plane z = a: (f)\ =0. In that case, for 

I z=a 

the propagator D{x,x'), ( ITOl) . which has to fulfill Dirichlet boundary conditions at 
z = a we take 

D{x, x') Dd{x, x') = D{x - x') - D{x - x') , (16) 
where x = (x, y, 2a — z, t) is the coordinate refiected on the plane z = a and 

/rl^h p-ik(x-x') 

is the free space propagator. Here m is the mass of the field but in the following 
we restrict ourselves to massless fields. 

In the next step we use the stationarity of the problem and turn to the Fourier 
transform in the time coordinate, X4 — > u. In this way the Casimir energy (|T5|l can 
be represented as 

1 /""^ 

E = — du trlnis:^ (18) 
2vr Jo 

with (cf. (ED) 

KUz,z') = DnAfi^)Ji^')) (19) 



5 



and 

Dd,M, x') = D^{x - x') -D^{x- x'). (20) 
as well as the free space propagator after Fourier transform in the time direction, 

/d^k Q-ik(x-x') 

In (fT8|) the trace is now only over the surface S and we used the symmetry of K^^ 
under u —u. We note that the function K^, ( fT9l) is even in u; — > —u. We take 
this into account in all integration over u below. 

Now let us consider a sphere of radius R centered in the origin so that the distance 
from its center to the plane is a. In terms of spherical coordinates (r, 6, ip) it is 
described by constant r and the angles provide a parameterization of S", 2; = {6, ip). 
A basis of functions in this space are the spherical harmonics | /,m) = Yi „^{6,(j)). 
Now we calculate the trace in f[T5]) in this basis, i.e., we represent the energy by 



^ = ^/ duj ili^rn^nKim^it^, (22) 







with 

Kim,i'm' = {l,m\ Dd{x,x') \l,m) 

and r^j. = r^./ = R. Now the trace is over the indices of the infinite dimensional 
matrix In Kim,i'm'- An equivalent representation is in terms of the determinant of 
this matrix, 

tTlnKim,l'm' =lndet{Kijn,i'm')- (23) 
Next we use the well known explicit expansion of the free space propagator, 

D^ix - x') = Y,yiA0cc,^.)di{r,r')Y;^,{9,,,ip,,), (24) 

with 

diir,r') = -^ji{ujr^)h\'\ury), (25) 

where ji{z) and h^i^\z) are the spherical Bessel functions to represent K in the form 
Ki,, = Si,, di{R, R) - {Im \D^{x- x') \ I'm) , (26) 
where because of the the azimuthal symmetry we could introduce 

Klm,l'm' = Sm,m'Ki,i'. (27) 
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The last step is to separate a distance independent contribution by writing 

InKi^v = In (Si,, di{R, R))+\n (^Si,, - ^^^^ ^^ {Im \ D^{i - x') \ I'm)^ . (28) 

The distance a enters this expression only trough x in the second term in the r.h.s. 
whereas the first is independent of a. In fact, it is this term which contains the 
surface divergences which are still present in the Casimir energy given by Eq. (I22l) . 
But dropping it we come to 

1 r°° 

E=— duj in^m In (5m' - ^i,v) , (29) 
where we introduced the notation 

= -J ^ D D\ (^"^ I ^'^(^ ~ I ^'"^) • (30) 
ai[n, n) 

The energy E, ( !29|) . is the ultraviolet finite, distance dependent part of the Casimir 
energy which be means of F = —^E is alone responsible for the force. 

The representation fl2^ is still to some extend symbolic. In fact, a precise formula 



is 



1 °° 

E = — diu V Indet(V-Ar), (31) 

Jo 



m=— oo 



where {6ij' — A;,/') is a infinite dimensional matrix with /,/' > |m|. 

Another way to represent the trace in an explicit formula is to expand the loga- 
rithm, 



^ -| OO OO OO oo 

E = ^ HvTi ^ E E ••• E (32) 

° 8=0 m=-(X) l=\m\h=\m\ ls=\m\ 

whereby the convergence of the sum over s is assumed. 

Let us mention that the representations fl3Tl) and fl32l) of the Casimir energy are 
exact. Their basic merit is that they are free of ultraviolet divergences. Hence they 
can be used also for direct numerical approximation. Examples of that are given in 
[6] and [7]. 

A further simplification appears because the matrix elements Ai^i can be written 
explicitly in terms of modified Bessel functions and 3j-symbols. In fact, the formula 



- V2V 2^ ^-^^ K,^.{uR) ^^^^ 

l+V / 

{-iY\2l" + l)K,.^.{2au:)i 



'I" I V \ f I" I I' 



i"=\i-i'\ 



/ V —m m 



holds. Its derivation is given in the Appendix A. 

Similar formulas as for the sphere can be written for a cylinder in front of a 
plane. For comparison we consider besides the (3+l)-dimensional also the (2+1)- 
dimensional circle in a front of a line on a plane. We denote the energies 

by E^l]_i-^ and -E|^2+i) correspondingly. 

We start from the formula ( |T8|) for the energy which for -E^^g+i) does not change, 

1 r°° 

Ef+i) = ^ / du trlnK^. (34) 

For (3+1) dimensions using the translational invariance along the axis of the cylinder 
it must be substituted by 

with 7 = ^Juj'^ + k1 where kz is the momentum parallel to the axis of the cylinder. 
Obviously the integration over can be done, 

1 /"°° 

^(^s+i) ^ ^ I du u tr\nK^ , (36) 

J 

so that the difference is just the additional factor ci;/2 in (l36i) as compared to ( l34l) . 
We note that E'^^^^^-^ is the energy per unit length of the cylinder. 

Next we have to consider the function K„, in two dimensions. A basis of functions 



on a circle is simply |m >= exp{imy9}/v2vr and the free space propagator has the 
decomposition 

+1) (3; _ a;') = ^ I ^) ^(2+1) I (37) 

m 

with 

d^^^^\r.r') = Ir^{ur^)K^{ur^). (38) 
In this way we obtain in place of 



Km,m' = 5m,m'dg+'\R, R) - {m \ oj^+'^i - x') \ m') , (39) 
where D^^'^^\x — x') is the same as (|2T1) in a dimension lower by one, 

/J2l. -ik{x-x') 

and X = {x, 2a — y) is the reflected coordinate in the (x,y)-plane. 

Separating again the distance independent contribution we represent the energy 
in the form 

^{3+1) ^ / ^ ^ln((5m,m/ - Am,ni') (41) 



with 

A^,m' = , I D^'^'^i^ - I ■ (42) 

In the (2+l)-dimensional case we have the same formula except for the factor u 
dropped, 



-^(2+1) 



— J du {6m,m' - ^m,mO • (43) 



Like in the case of a sphere, the matrix elements Am^m' can be calculated explic- 
itly (see Appendix A) and expressed in terms of Bessel functions, 

Am,m' = -flA rKm+m'C^UJa^ImiuJr) (44) 

a formula which coincides with the corresponding ones in [7] and [6] . 



4 Cylinder at small separation from a plane 

In this section we consider a cylinder at small separation from a plane. We start with 
Dirichlet boundary conditions. Here we consider first the case of (3+1) dimensions, 
i.e., the usual space. Then we reduce the dimension by one, i.e., we consider the 
(2+l)-dimensional case which is a circle in front of a line on a plane. Finally, we 
consider the case of Neumann boundary conditions. Since in the cylindrical case the 
polarizations of the electromagnetic field separate into TM modes fulfilling Dirichlet 
boundary conditions and TE modes fulfilling Neumann boundary conditions this 
gives at once the result for the electromagnetic field. 



4.1 Dirichlet case 

Starting point is the formula fHTj) . We rescale uj —>■ uj/r so that the energy is now 
given by 

1 

^{3+1) = / ^ tl In {Sm,m' - Am,ni'{^^)) , (45) 

where Am,m' is now given by 

(l^) = I N Km+m'{'2'Uj{l + e)) /m'(i^) (46) 

with 

where L = a — R is the separation between the cylinder and the plane and we are 
interested in the small distance behavior, i.e., in e — * 0. 
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To start with, we expand the logarithm in Eq. (l45|) . Below this turns out to be 
justified because in the orders in e we are interested in this series converges. The 
energy is now given by 

s=l m ni ris 

with 

= ^m,m+nx{^)^m+nx,m+n2{^) ■ ■ ■ ^m+ns-i,m+ns{^) ^m+na,m{^) ■ (49) 

In ( l48l) the summations over m and the n^'s run over the integers. 

Now the small distance behavior of the energy comes from the region in the 
integration resp. sums of large all, cj, m and the n^'s. We use the uniform asymptotic 
expansion of the Bessel functions and substitute the sums by integrals. Further 
we use the symmetry under m —>■ —m and obtain in the sense of an asymptotic 
expansion 

^ -j^ poo poo poo POO 

where for the Am^m'i^) in -M. still Eq. (l46l) holds, however, with the Bessel functions 
replaced by their well known uniform asymptotic expansions which are displayed in 
the Appendix B. From the expressions in the exponential factors in 

Am,m' ~ e"^ , (51) 

which combine into 

, /2a;(l + e)\ fuj\ , f uj \ 

7] = [m + m j 1] [ \ — m 1] [ — )— m rj [ — ) (52) 

\ m + m' J \mJ \m' / 

(see Eqs.( l46l) and ( IB.ll) ). it follows that for small e the dominating contribution into 
( l50l) comes from ti;~ m~i, nj~^. Hence we substitute 

tVl - r2 tT 

uj = , m = — , 

e e 

Hi Hi d—, (53) 

where the factor y/^i was introduced for later convenience. With these substitutions 
the energy becomes 

_ -3 °° 1 poo r, nl i>00 POO / A+\ i 
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Now we change the notations for the Am,m'{^^), 

indicating the dependence on the variables n and n'. Using the asymptotic expan- 
sions of the Bessel functions entering by means of Eq. fHB]) . we obtain the expansion 



<n' = ^^^ e-'^ (l + v^<i'^(t,^) + e<Ut,r) + ...) (56) 

with 

ri"^ = 2t+{n- n'f. (57) 

Details including the explicit form of the functions a^^^^^(t,r) and (I'nn'i^i'^) 
shown in Appendix B. 

Next we insert (|56|) for the A^.m^"^) i^to (H9l) and after a reexpansion we 

obtain 



\0<i<j<s i=0 / 

with 

s 

?7i = (nj - nj+i)^ (59) 

Here, in order to include all contributions into the sum signs, we have to put uq = 
Us+i = formally. 

The above formulas allow us to to represent the energy E, Eq. flM|) . in the form 

(60) 

with 

= e-^^ (l + Ve Y.a^^;X.{t,r) (61) 
+^ ( E . + E ^U. it,r))+..). 

\0<i<j<s i=0 / ) 

In this representation it is possible to carry out the integrations over r and over t 
and we obtain 

\ \ dm p dus r r(5/2) 



i=0 



^(3+1) 27rL2Vi:^S + l./_^^/^••■7_^^/i 



47r(2(s + l))5/2 
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b{ni, ni+i, nj,nj+i) + ^ d{ni, n^+i) + . . . | + . . . V (62) 



with 



b{ni,ni+i,nj,nj+i) = 

dt tfe-2(-+i)* 



dr o!^!X,it.r)c^j%^^.r) (63) 



and 



t v47r 7n 



d(n,,n,+i) = / y y dT a^^l^^^^{t,T) . (64) 



The quite lengthy exphcit formulas for h and d are shown at the end of Appendix 
B, Eqs. f lB.14p and fIB.lSp . In Eq. fl62l) we dropped the contribution proportional to 
y/e. It can be seen that it contains only odd powers of the n^'s and vanishes under 
the integration over the n^'s. 

The next step is to perform the integrations over the n^'s in Eq.( l62l) . In the 
leading order in e this is quite simple. We rewrite the quadratic form r^i, Eq. (!59!l . 
by noticing that obviously the formula 

H 7 n^-i H (65) 

s — 1 \ s J s 

holds. It immediately delivers 

drii f°° dn 



00 



_ e""^ = , (66) 



and we obtain the energy for a cylinder in front of a plane in leading order in small 
separation, 

pcyl PFT ^ ^ ^ r(5/2) 2-^/^ 



s=0 



-1 fR 3C(4) _ TT^ 1 /i? ^^^^ 



27rL2 V L 32v^ 1920v^^2V L ' 

which coincides with the known result from applying the proximity force theorem. 
Note that in order to restore dimensions it has to be multiplied with h and c. Also 
note that this is the energy per unit length of the cylinder. 
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Now we turn to the remaining terms in Eg. (162!) . Of course, all integrations there 
over the nj's are Gaussian. But the combinatorics is somewhat involved so that we 
delegate these calculations into Appendix C. The result is enjoyable simple, 



drii 



-oo ^' 



I -7= \ ^ b{ni,ni+i,nj,nj+i) +^d{ 

J-00 \0<i<j<s i=0 



384^2 (1 + s) 



The remaining sum over s in (!62l) is the same as in the leading order and finally for 
small distance L the energy becomes 



''(3+i)-27rL2VL 32v^ 1 ^QR \[rJ 



2^ 



(69) 



a formula, which represents the first correction to the proximity force theorem. 

Using the method shown above it is easy to repeat the calculation for (2+1) 
dimensions. The first change is in the substitution (1531) where it is now meaningful 
not to use the variable r but to define 

t sin 09 t cos if 

oj = , m = , 

e e 

Ui rii (70) 
with Lf G [0, |]. Next we mention the changed formula (I6OI) 



as 



(2+1)" 2lV I ^7 + 170 T Jo ^ J-ooT^'" J-ooV^ 

(71) 

where besides the prefactor only the power of t changed and the integration over 
T is changed for that over cp. The remaining steps go in complete parallel to the 
(3+l)-dimensional case and the result is 



4.2 Neumann Case 

Let us consider the changes which come in for Neumann boundary conditions first 
in the general formulas of section 2. Basically it is sufficient to change in formula ([6]) 
the delta function for its normal (to the surface S) derivative, H dn^6{x — f{z)). 
As a consequence, the definition of K{z,z'), ([9]), changes for 

K{z,z') = dnA^,D{fix)J{x')), (73) 
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where the normal derivatives must be taken before putting the argument onto the 
surface. The next changes is in formula (fT6!) because on the plane at 2 = a we have 
to take Neumann conditions too. It amounts simply in a different sign, 

D{x, x) Dn{x, x) = D{x — x) + D{x — x) . (74) 

As next changes we mention formula (l38l) where the Bessel functions must be sub- 
stituted by their derivatives, 

dt'-'\ry) = r^{ur^)K'^{ur^). (75) 

and equation (14T|) which looks now 

^(I+l) ^ ^ j ^ ^ln{6m,m' +Am,m') (76) 

with ^ 

Am,m' = ---—-rK^+m'{2uja)r^{ujR) (77) 

Now we trace the changes in the asymptotic expansion for small separations. First we 
consider the asymptotic expansion of Am,m', fTTTj) . using the asymptotic expansions 
of the derivatives of the Bessel functions, 

I'Jmz) \_^n^Ul + ^..,(.)) U^-_M^,\ (78) 



K'^imz) J ^J2m z \ m 

with Vi{t) = (— 9t + 7t'^)/24. The first observation is that the minus sign in the 
expansion of K'^{mz) just compensates the sign change in (17^ . The second remark 
is that in the expansion of the factors in front of the exponential, (1B.6P in the 
Dirichlet case, just n and n' change their places which does not affect the expression 
under the trace in fl76l) . The final remark is that the only change comes from the 
Debye polynomial Vi{t) which appears in two places instead of Ui{t). The further 
calculation runs then in exactly the same way as in the Dirichlet case, but the result 
is different, namely for the formula which corresponds to (jHSD we obtain 



dni 
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I ^ / ^ 6(ni,ni+i,nj,nj+i) +^c?( 

J -co \0<i<j<s 1=0 



(79) 



384^2(1 + 5)3 12v^(s + l) ' 

(80) 

i.e., we get a contribution which is the same as in the Dirichlet case and an additional 
contribution. Inserted into the energy this gives 



/^3C(4) f ^1_8C(2)\ 
■L^^ L 32v^\ V36 9C(4)y' 



2^ 
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a formula, which represents the first correction to the proximity force theorem for 
Neumann boundary conditions. 

We note that the mentioned compensation of the signs ensures that the leading 
order contribution, i.e., the proximity force theorem, gives the same answer for 
both types of boundary conditions as expected. Less obvious is the behavior of the 
Neumann contribution in what it has in addition to the Dirichlet case (which equals 
the first term in the r.h.s. of Eq. fl79|) ). This is the second term in the r.h.s. of Eq. 
f l79|) . which results in a much less convergent sum over s. In the (3+1) dimensional 
case the sum is still convergent, but in the (2+1) dimensional case not. Carrying 
out the corresponding calculation we obtain 



32^2(5 + 1)2 4^2 

in place of the r.h.s. of Eq.( 179i) . Now the sum over s is logarithmic divergent. In fact 
this implies that the correction to the proximity force theorem in (2+1) dimensions 
for Neumann boundary conditions is proportional to ^ln(L/i?). 

Finally we collect the energies for Dirichlet and Neumann boundary conditions 
together to the first correction to the proximity force theorem for the electromagnetic 
fields. 



T^cyl,electromgn. 
-^(3+1) 



-1 [RsgA) j 

27rL2 V L 16^2 \ 
1 

~Z2 




^3) 



5 Sphere and cylinder at large separation from a 
plane 

The limiting case of large separation between a sphere and a cylinder and a plane 
is known and we display it here for completeness. For a sphere the starting point is 
Eq. (!32!) with Ai^i given by Eq. (l33!) . The main contribution comes from small ut and 
small orbital momenta. We substitute 

. - ^ (84) 

and expand Ai i/ for small R/a. simply by expanding the two Bessel functions in 
fl33|) which depend on r. The leading order is 

R 

Ai,u = Si,o6i',o ^ e"" + . . . . (85) 
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Inserted into E, (l32l) . this gives immediately the large distance limit 




where we included the first correction too (note a = L + R). This result agrees with 

The corresponding procedure for a cylinder is similar. We start from Eq.( 
or from Eq. fHS]) in dependence on the dimension and with Am,m' given by Eq. 
Again, the dominating contribution comes from small momenta. We do the substi- 
tution (18^ and expand Am,m'- Here, however, a peculiarity shows up because of the 
logarithmic behavior of the Bessel function Kmi^^-^) in the denominator for m = 0. 
The expansion is 

^^,^' = ^^^^#^^ + 01 (-) 1 , (87) 




-7 - \n{ujR/Aa) 

where 7 is Euler's constant. Inserting into fl32l) and integrating over u we get 



(2+1) 8a \n{R/2a) \ \n{R/2a 
and 




7rVl2 

(3+1) - ln{R/2a) (ln(i?/2a))2 ' ^ \ \ln{R/2a) 



= — — 1 1 - :j,r... + o 




9) 



In the (3+l)-dimensional case the next-to-leading order is zero. 

We note that the logarithmic behavior of the large distance limit is related to 
the corresponding behavior of the Greens function in two dimensions. Because the 
logarithmic behavior is present in E^^l^-^^^ i.e., for a circle on a plane, too it seems 
unlikely that it is related to the length of the cylinder as discussed for the three 
dimensional case in [7]. 



6 Conclusions 

In the foregoing sections we showed the derivation of the expression for the Casimir 
energy in terms of the determinant of the projection K{z, z') of the free space prop- 
agator on the boundary surface S which was known from [H] and [TU] and indepen- 
dently rederived and applied to the calculation at large distance in [6] and [7] which 
can be considered as a significant progress because it opened the way for direct nu- 
merical calculation of the Casimir force for non simple geometries. In the present 
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paper we showed that also the short distance behavior can be calculated by this 
method and derived the first corrections to the force proximity theorem including 
the electromagnetic case. In (3+1) dimensions the corrections are power like, in 
(2+1) dimensions the corrections contain a logarithm for Neumann conditions. 

In general, corrections to the limiting behavior for both, large and small dis- 
tance, can be calculated as expansion in the corresponding small parameter. These 
expansions are asymptotic ones. For the small distance behavior this is clear from 
section 4 where we substituted the sums over the orbital momenta of the photon 
propagator by the corresponding integrals. Obviously, exponentially small contri- 
butions are neglected in this way. For the long distance expansion this follows from 
the denominators in the corresponding matrix elements, (l33l) and (jH]). Consider the 
energy as a function of the expansion parameter 5 = r/2a (after the rescaling (!84|) ) in 
the complex plane of 5. When moving 5 around the origin the pole of the integrand 
crosses at some time the path of the integration over uj delivering an additional 
contribution. Hence, the energy as a function of 5 has a cut starting from (5 = 0. 

The new method for the Casimir force for non simple geometry which so far 
gave the possibility to calculate the large and the small distance behavior, has the 
potential in it for the calculation of the Casimir force at all distances with desired 
precision. Also, it obviously can be generalized to much more boundary conditions 
than Dirichlet or Neumann ones. In the present paper we considered primarily 
Dirichlet conditions in order to keep the representation as simple as possible. The 
generalizations to Neumann are obvious, one has to change the sign in (fT6!) and to 
add a derivative to H{x, z), Eq.Q as done, for example, in [7j. In the path integral 
approach this was considered from a more general point of view in |llj . 

The short distance expansion in section 4 sheds light also on the quasi classic, op- 
tical path and similar approximations. Indeed, the dominating contributions come 
solely from large momentum, however in an expansion over the number of reflec- 
tions, like the expansion of the logarithm in section 4 shows, all terms contribute. 
Therefore, in order to get some approximation, say in the first order of the small 
parameter, all reflections must be included. 
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Appendix A 

In this appendix we show the calculation of the matrix elements ( H2l) and ( l30i) in 
section 3. We start with the cylindrical case. Using (HOj) we represent the matrix 
element in ( 1421) in the form 

W^^m' = {m\D^J+'\x-x') \m') 

Ml. p2iak2 

fmCk) fm (A.l) 



2 



with 



(27r)2 a;2 + k 







Here we used x = {xi,2a — X2) and defined k = (fci, — ^2) in order to rewrite 
k{x — x') = 2ak2 + kx — kx'. Using polar coordinates in the fc-plane, ki = ksimpi^, 
k2 = kcos(fk, and an integral representation 

6"^™^+^""'"^ = Jm{z) (A.3) 

271 

of the Bessel function Jm{z) for the functions fm{k) entering ( lA.ip we obtain 

/J2r, 2iak2 
^ -^T7^ '^mikr) Jm'ikr) e-^('"+'" (A.4) 
ivr uj^ + k^ 







The formula (lA.Sp can also be used for the integration over the angle (pk, 



dk k 

Wm,m' = I ^ , ,p Jmikr) Jm'ikr) Jm+m'i^ka). (A.5) 







CU2 + F 



The last integration can be done by extending it over the whole /c-axis and clos- 
ing the contour picking up the pole in k = iuo. With the relations Jm{z) = 

f(^H^\z)+H^\z)'), Hi^\-z) = i-ir+^H!^\z), M-z) = i-irM^) one 

obtains 

Wm,m' = — Jmiii^r) Jm'(it^r) (2it^a) , (A.6) 

or in terms of the modified Bessel functions, 

Wm ,m' 

which after division by dm'^^\r,r), (138!) . gives (jHj). 

Now we turn to the spherical case and consider the matrix elements in (l30l) . 
Using (I2TI) we represent them in the form 

Di^i> = {Im \ D^{x - x) \ I'm) 
18 



with k = {ki, k2, —k-s) and 



2n 



Jo 

= A7i{-tyjiikr)YU9k,Vk) , (A.9) 
whereby in the last hne the decomposition of a plane wave into spherical waves, 

e^'^ = 47r^2'jKA;r)l^:„(^„^,)y^„(4,^,) , (A.IO) 



Im 



was used. Here, ji{z) are the spherical Bessel functions. Writing the integration over 
k in spherical coordinates and using again ( lA.lOl) an integral over a product of three 
hypergeometric functions appears which can be written in terms of 3j-symbols, 



/' 

^0 



r2n 

d9sm6 I dip Yirn{0,(p)Yii„,'{9,(p)Yi>i„,„{e,(p) 



(2/ + l)(2/' + l)(2F + l) [ I V I" \ [ I I' I" 



An \ J \ m m' m" 

and we arrive at 



(A.ll^ 



Di,, = -v/(2/ + l)(2Z' + l) / -^^j,{kr)Mkr) (A.12) 

E >'--'"(2'"+ibH2.^-)(;;' ^ i 



i"=\i~i'\ 



Using ji{z) = ^{h\^\z) + hf'\z))^ the integration can be extended to the whole 
/c-axis and closed picking up the pole 'mk = iuj. The result can be inserted into (130|) 
and we obtain 



\v = V(2/ + l)(2/- + l) /;;|'''''^^ (A. 13) 

i+i' 



J2 i'^''-'"{2l" + l)h'i^}{2iau 
i"=\i-i'\ 



I" I V \ f I" I v 



/ V — m m 



which coincides with the corresponding formula in [6J. Passing now to the modified 
Bessel functions and dividing by di{r,r), fl25l) . we obtain just formula (|33|1 . 
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Appendix B 

In this appendix we calculate the uniform asymptotic expansion of Am,Tn'{^), Eq.( l46i) . 
using (IB.ip and obtain A'^^,, Eq. (!56|) . and the coefficients all^^\t, r) and cLn\'i^i 
therein. 

We remind the uniform asymptotic expansion of the Bessel functions, 



Kmimz) \~J2^ (1 + ^2)1/4 



+ (B.ll 

m 



wit h = VTT^ + ln(z/(l + VTT^)) and Ui{t) = (3t - 5t=^)/24 with t = 
1 / ^/l + z^, the first of the Debye polynomials. Since we intend to calculate the first 
next-to- leading order only the restriction to ui{t) is sufficient. 

We start with the exponential factor and do the substitution (|53l) there. Up to 
order e we get 

fj = 2t + {n~ n' f + JjT (2t - [n - n'f) [n + n') (B.2) 



+e ( -r H ^ — [n + n) 12^^ ((n-n) + 6(n - n 

The expansion of the exponential is 

= (l + v'tln^ + 6 , + 0)e^)) (B.3) 

with 



V^nln^ = T (2t -in- n'r) {n + n'). (B.4) 



and 



ryW, = _i + (1 - T^){n + n')' - ^J^l^ " n'f {In'' + lOW + 7n''). (B.5) 

Next we consider the factors in front of the exponential in A^^,. From (146|) and 
( IB. II) they collect into 

1 /m 1 + 1 ^^ 



V 27r \m' m + m' J \1 + Z2I + z^J 

with zi = —, Z2 = and 2:1 = Ml±£l_ Here, again, we do the substitution fl351) and 
the expansion for small e gives 
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with 



and 



.(1) 



+ 1 



Qn^ - Ann - lOra'^ + {-lln^ + 2nn + 2Qn^) ) . (B.9) 



(B.8) 



/2n 2 



We note hat these contributions are not symmetric under n n'. 

Finally we have to consider the contributions from the Debye polynomials. To 
the given order we note 



2 _|_ _j_ 
D = — 

1 _ Ui{t2) I 

^ m' ^ ■ ■ ■ 



2 Mh) ^ 
m + m' 



(B.IO) 



with ti = a/1 + zf. Doing here the substitution (1531) we obtain 

D-^l + e + ... 

with 



16t 



(B.ll) 



(B.12) 



Taking flB.3l) . flB.7p and fIB.llI) together and reexpanding again we obtain finally for 
a';lff{t,T) and a«,(t,r) in Eq.® 



{n - 3n' + 2{n + n') [n^ - 2n'n + n'^ - 2t)) r 

{n - d,n') {n + n') {n^ - 2n'n + n'^ - 2t) 3 - 5r 



2t 



16t 



+ 



+ 



8t 



I2t 



(6 - llr^) + 2n' (r^ - 2) n + 4t (r^ - l) 
+ n'^ (29r2 - lO) 
"-7(3r2-l)n^ + 4n' (3x^-1) 
+ 6 (^(3r2 - 1) n'^ + 2t (r^ - l) j ^ 
+ W (^(3r2 - 1) n'^ + 6t (r^ - l) 
+ 6(n + n )^ (^n^ - 2nn + n'^ - 2t j ^ 
+ 12tV2 + I2n'h [t^ - 1) - 7n'' (3r2 - l) 1 . (B.13) 
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At the end of this appendix we display the exphcit forms of the functions b, (l63l) . 
and d, (164|) . introduced in Eq. (!62|) : 



b(x, y, z, w) = ^ l"4(s + 1) (4(s + l)x^ - 4(s + l)yx'^ 

+2 (-2(s + l)y^ + s-2)x + y (4(s + l)y^ - 6(s + 1) - 6)) 

+4(s + 1) (-4(s + l)x^ + 4(s + l)yx^ + (4(s + l)y^ - 2{s + 1) + 6) a; 

+y (-4(s + l)y^ + Q{s + 1) + 6)) zw^ + (-4(s + 1) {2{s + 1) {2z^ + 3) + 6) 

+4(s + 1)1/ (2(s + 1) (2^2 + 3) + g) ^2 ^ ^4 ^2|/2 _ 1) (2^2 + 3) + i)2 

+ 24 + ^2 ^ ^ ^) ^ gQ^ ^ ^ ^ (^_4 (^2y2 _ 3) (2^2 + 3) + 1)2 

-24 {y^ -z^ -3){s + l) + 60)) w + z (8(s + 1) {2sz^ + 2z^ + s - 2) x^ 
-8{s + l)y {2sz^ + 2z'^ + s-2)x'^ + (-4 {2y'^ - l) {2z'^ + l) (s + 1)^ 
+24 _ ^2 _ ^ ^ gg^ ^ ^ ^ ^4 (^2|/2 _ 3) {2z^ + 1) (s + 1)2 

-24 (y' + Z-l) (s + l) + 60)) 



(B.14) 



1 



(s + 1)V 



d(x,y) = ■= 

^ ^ 384v^(s + 1)7/2 

-16(s + ifyx^ - 8(s + 1) ((s + l)y^ -3 + 2)2 

+32(s + ify {y^ - l) x^ - 2 (4(s + l)y - 8(s 

+ (4 - 7s)s - 4) + 4y (-4(s + 1) V + 8(s + 1) V - s(5s + 16) + 4) x 

+8(s + 1) V - 24(s + l)(s + 2)y^ + (40 - 2(s - A)s)y^ + 4(s - l)s + t] . (B.15) 



+ l){s + 4)y' 



Appendix C 

In this appendix we calculate the integrals over the n^'s in Eq.f l62p . As mentioned in 
the text all these integrations are Gaussian and we are left with the corresponding 
combinatorics. 

We note that a number of these integrations can be carried out independently 
of the functions 6(nj, nj+i, nj, n^+i) and d{ni,ni+i) entering (1^ . In the following 
we divide the whole expression into a number of pieces and in each we first perform 
these integrations. The remaining ones, one to four in number, depend on the details 
of these functions. They have been performed machined. 

We start with the part involving d{ni,ni+i), (|6^ . in fl62l) . It is easier because it 
contains only one sum. First we separate the first and the last contributions to this 
sum, 

s s—l 

d{ni, rii+i) = d{0, rii) + d{ns, 0) + ^ d{ni, n^+i). (C.l) 

4 = 4 = 1 

Renaming the nj's by numbering them in the reverse order in the first term in the 
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r.h.s. we have 



S, ^ I ^ . . . r ^ (d(0, m) + d{n,, 0)) e-"^ 



— oo 



= / ^.../ ^(rf(0,n,) + dKO)) e-^\ (C.2) 

Now we apply formula ( l65l) . However, in order to avoid confusion with the notations 
we rewrite it, 



ml + Y^ {rrii - rrii+if + ml =j (^nii - ^^2^ + ^ (^^2 - ^mg^ + . . . 



• • • H rriu H (C.3) 

u — 1 \ u J u 

With M = s, mi = rii, m„_i = n^.i the first (s — 1) integrations can be carried 
out, one by one. The first is that over ni, which must be shifted, ni — > rii + 

and gives just y^v^- The second integrations needs the shift ?22 — > ^2 + f'^s and 

gives -y/l and so on. The shifts of the variables which do not enter the functions 

d (and, below, h) are obvious and we do not mention them in the following. Finally, 
renaming = x we obtain 

^ f°° dr / s -I- 1 \ 

Si = ^ UX ^^^^^ ^ ^^^^ _ (Q_4) 

The last integration was carried out machined (like the similar ones below) and the 
result reads 

S - 1 3 I 37 1__ 

' 12y2(s + l)2 8v^(s + l)3 64v^(s + l)4 4v^(s + l)5' 

Next we calculate the contribution from the sum over i in (IC.1|) . To this end we 
represent rji in the form 

Vi = ^1 + (^1 - ^2)^ + . . . + (ni_i - riiY + 

2 I / \2 2 

-n^ + [rii - rii+i) - n^^^^ 

+ (rii+i - ni^2? + ■ ■ ■ + (C.6) 

and apply fIC.Sp to ni, rti with u = z and mi = rii, m„ = rij and to rij+i, 
Hs with u = s — i and mi = n^, m^j = nj+i so that it takes the form 

2 / 1 \^ if i-l " 

2/ s — 1 \^ s — i ( s — i — \ 



+ T + . . . H ^ — ^i+i 

1\ 2 / s — ? — 1\ s — « 
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With this formula the integrations over ni, and over nj+i, Us can be done 

and renaming rii = x and nj+i = y we obtain 

drii f°° dus ^ ^ 



S2 ^ / / ^ (C. 



i + 1 / i \ s + 1 



2 



■ exp — X - — — y - -y 

The last two integrations and the sum over i give 

^ 16^2(5 + 1)2 384v^(s + l)3 192^2(5 + 1)4 472(5 + 1)5' 



Taken together the expressions simplify to some extend 



1 17 1 

Si + 82 = ^ ^ + ^ . (C.IO) 

48v^(5+l)2 384^2(5 + 1)3 24^2(5 + 1)4 

Now we calculate the contributions of the double sum in (!62l) . Again we separate 
the first and the last terms, 



^-00 V -'-00 V 0<i<i<s 

dn\ /"°° dUs ( , „^ r \ 

\ — 6(0, ni, 72^,0) + > 6(0,ni,nj,nj+i) 



0<j<s 

+ ^ 6(ni,ni+i,ns,0) + ^ 6(ni, n^+i, n^, n^+i) j 

0<j<s 0<i<i<s 

= 'S's^i + + 'S's^a + S'3_4. (C.ll) 

For ^3 1 we rewrite 771 in the form 

r]\ = nl + {n2 - riif + {n2 - rii - {n3 - ni)f + . . . 

+ (ris-i - rii - {us - ni)f + (ra^ - rii)^ 
-(^.-m)' + n2 (C.12) 

and apply flC.31) with m = 5 — 1, mi = n2 — rii, = — rii, 

2 / 1 , ,V 1 s-2, ' ' 

Vi = J {^2 - ni - -{ns - ni) j + ■ ■ ■ + j—^ {Us^i - Ui - -—j{ns - Ui) 

ni--n)l +'^nl (C.13) 



5 — 1 V 5 / 5 
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After integration over n2,.--, "^s-i and renaming ni = x and Ug = y we have 
83,1 = I J / / ^o(0,x,?/,0) 



The last two integrations give 



^3 1 = 7^ ^ • (C.15) 

For ^3^2 we rewrite ri\ in the form 

r\\ = nl + (ni - ^2)^ + . . . + - n^)^ + ra^ 
-ra^^ + (ui - Ui+iY - (ui+i - n,f 
+ {Ui+i - Tisf + (n^+i - - (ni+2 - ns)f + ... 

+ (ras-2 -ris- {us-i - ns)f + {us-i - risf + nl (C.16) 



and apply fIC.SP to rii, rij with u = i, mi = rii, = and to rij+i — ra^, 

Us-i — Us with u = s — i — 1 and mi = nj+i — Ug, rriu = Ug-i — Ug so that rji 
takes the form 

+Y I n^-i - - 2 (^^-2 - ^s) I + • • • 

s — i — 1 f s — i — 2 

H : ni+2 - ris : rl^^i+i - ris) 

s — I — 2, \ s — I — 1 

•,1/ • \2 / ■,-|\2 

1 + 1 I \ S I l + 1 



rii - — -Tij+i + 7 : , Ui+i n. 



' I V ' ^+1 V (s-z-l)(z + l) V ' s 

^5 + 1 2 

H n, . 

s 

After integration over rii,..., r;,j_i and nj+i,..., n^-i and renaming jt-i = x, nj+i = ?/ 
and ris = z we have 

1 P ^^a; p rfy p ^ 

•53,2 = ^ . = / ^/ ^/ ^6(^,Z/,z,0) (C.18) 
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The last three integrations and the sum give 

In the same way we calculate 5*3^3 . Here we rewrite rji in the form 

Vi = + (n2 - riif + (722 - rii - (^3 - ni)f . . . 

+ {rij^i -rii - {rij - ni)f + {nj - riif + 
-{uj -nif + {rij - rij+if - n]^-^ 

+n]^^ + K+i - n,+2))' + ... + nl (C.20) 

and apply fIC.Sp to ^2, nj_i with u = j — 1, mi = ^2 — "^i, m„ = Uj^i — rii 
and to rij+i, with u = s — j and mi = ng, = 'n.j+i so that r/i takes the 

form 

= Y I -^1 - 2('^3 -^1) I + ■ ■ ■ + \nj_i - ni - -—j{nj - ni) 

2/ 1 V s-j / s-j-1 

IV 2 J s-j -1 \ s-j 

S+1 2 S ( S-J+l 

H — — rij rii 

s — j -\- 1. f s — j \ ^ 



S-J V 

After integration over n2,---, nj^i and over nj+2v) and renaming ni = x, rij = y 
and Hg = z we have 

1 P f^a; r dy r dz , 

-^3,3 = E^7=^T^=Y / ^ / ^ / ^KO,x,i/,z C.22 

-PI-—- -(,-,H-l)U-l) i^' 

+ (C.23) 



The last three integrations and the sum give 

^3 3 = IZI . (C.24) 

8^2(5 + 1)4 ^ ^ 

The last step is 5*3 4, which in fact contains the double sum. Here we represent 
rji in the form 

rii = nl + {rii - 722)^ + . . . + [rii^i - riif + 
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-n- + {rii - Ui+if 

+ (^i+2 - ni+iY + ini+2 - iii+i - - n^+i))^ + . . . 

+ {nj+i - nj+2)f + . . . + n^. (C.25) 

We apply (1C.3P three times, first, with u = i and mi = rii, = n^, second, with 

u = j — i — 1 and mi = nj+2 — rij+i, niu = rij — n^+i, and, third, with u = s — j 
and mi = rig, m^ = rij+i. So we write rji in the form 



1 V z 
^1 no + . . . H 

2 / i-1 



i - 1 
—ni 



+ - ni+2 - ni+i - -[rii+s - rij+i) + . . . + 



+ 



j - ^ - 1 
j-t-2 



?2j_i - nj+i - -. -{rij - ni+2 



J 



+ T \ns- -ns-i + . . . + 



s-j-1 



i + 1 



rii - — pY^i+i ) + 



S + 1 2 S - j + 1 

j(s - J + 1) s- 2 



s — j — 1 

nj+2 — '^i+i 

s- 3 



i + l 
rii+i — % 

g — J \ 



All integration except for four, whose variables we rename for rii = x, rii+i = y, 
rij = z and rij+i = ty, can be done. 



3,4 



E 



1 



0<i<j<; 



rfx /""^ c/y /""^ dz dw 

-1= b{^,y,z,w} 



■ exp 



J — oo 
Z + 1 



— oo y ^ J — oo 
2 



i + V 



i + l 

y —z 



s + l 



{s -] + !)]' 
The last integrations give 



s-3 + l 
s-3 



w 



2 H 



(C.27) 



(2z + l)(4j-3s-l) 



(C.28) 
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and after the summations over i and j 



48v^(s + l)5 ^ ^ 



Together we get for ^3, flOTTD . 



A8V2{s + 1)4 ■ 
Collecting together all parts we obtain finally 



S, = -J^^- (C.30) 



^1 + ^2 + ^3 = 7^ , (C.31) 

384^2(5 + 1)3 ' ^ ^ 

which is the result quoted in Eq. fl68|) . 
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